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Abstract - We discuss zeta functions, and traces of the associated weighted evolution 
^ ■ operators for intermittent Hamiltonian systems in general and for the Sinai billiard 

I in particular. The intermittency of this billiard is utilized so that the zeta functions 

■ may be approximately expressed in terms of the probability distribution of laminar 

I lenghts. In particular we study a one-parameter family of weights. Depending on 

• the parameter the trace can be dominated by branch cuts in the zeta function or by 

I isolated zeros. In the former case the time dependence of the trace is dominated by 

0^ ' a powerlaw and in the latter case by an exponential. A phase transition occurs when 

the leading zero collides with a branch cut. The family considered is relevant for the 
^1^. calculation of resonance spectra, semiclassical spectra and topological entropy. 

o 

1 Introduction 

, The space of Hamiltonian systems is framed by two extreme cases - integrable systems and 

^ ' strongly chaotic systems, the latter in the sense of Axiom-A. Both cases may be successfully 

solved, but the methods of solution look entirely different. For integrable systems solving means 
finding the constants of motion. In the Axiom-A case one rather aims at a statistic description 
of the motion - one calculates different kinds of entropies, fractal dimensions, resonance spectra 
etc. To this end cycle expansions of zeta functions have shown very successful |, ||. So far, 
the success of this approach has been limited to the Axiom-A case. The reason is that Axiom-A 
ensures nice analytic features of the zeta functions involved [§[ ||, , whereas different kinds of 
singularities emerges when departing from Axiom-A If the system is not Axiom-A, but still 
ergodic, it is generally intermittent (there are cycles with arbitrary small Lyapunov exponents) 
and lacks a simple symbolic dynamics. This is e.g. the case for most bound ergodic systems 
It is of great interest to extend the study of zeta functions to this case. The resulting theory 
would be an analogue of perturbation theory for almost integrable systems. 

Such a theory would be interesting not only for ergodic systems. The motion in a chaotic 
sea of a generic mixed system is in general intermittent. This is because the trajectory may be 
trapped inside cantori surrounding the stable islands, where it exhibits quasi-integrable motion. 

The cycle expansion for a system with a finite subshift symbolic dynamics consists of two parts 
jl] , a fundamental part, giving the gross structure of the spectrum, and curvature corrections. 

When the symbolic dynamics is an infinite subshift, there is no similar division. One goal 
of this paper is to establish a fundamental part for this case. One can think of two ways out 
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of the dilemma. One can define succesive finite subshifts, approaching the infinite subshift, by 
utihzing a pruning front |jT^, . This mean that one persists on a periodic orbit description. We 
will attack the problem from the opposite end and completely abandon the periodic orbits. We 
will thus make extensive use of our systems being ergodic and Hamiltonian, so that the invariant 
density is known a priori. One of the morals to be taught is that there is no reason to find a lot 
of periodic orbits in order to compute something you already know in advance. The idea to use 
probabilities rather than periodic orbits has been used in a similar way in ref. . 

In section 2 we develope most of the theoretical apparatus needed. The exposition will be 
rather brief as most of the material can be found elsewhere. In section 3, we apply these ideas 
to the Sinai billiard and formulate an approximate zeta function. In section 4 we study this zeta 
function numerically and compare with periodic orbit theory. 



2 Theory 

Much of the theoretical work on these matters are centred around the evolution operator. It 
describes the evolution of a phase space density <I>(a;) 

£mx)= J wix,t)5{x-f{y)My)dy . (1) 

The phase space point x is taken by the flow to f*{x) during time t. w{x,t) is a weight 
associated with a trajectory starting at x and evolved during time t. It is multiplicative along 
the flow, that is w{x, ti + ~ w{x, ti)w{f*^ {x), 12)- 

We emphasize that i is a continous variable as we are studying flows and not maps. 



2.1 Trace formulas and zeta functions 

Our main concern in this paper is to compute the trace of the evolution operator, that is, the 
sum of its eigenvalues. The material in this subsection may be found e.g. in refs. [|l|, |l^ 
The trace may be written as a sum over the periodic orbits in the system 



trC}^^ I ■w{x,t)5{x - f{x))dx 
5{t-nTp) 
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Z.^pZ.^P|det(l-M")| ' 

p n— 1 ^ 

where n is the number of repetitions of primitive orbit p, having period Tp, and Mp is the Jacobian 
of the Poincare map. Wp is the weight associated with cycle p. 

The trace is written as the Fourier transform of the logarithmic derivative of a zeta function 

_ 1 /"°° Akt ^'wi^) . 



We restrict ourselves to systems with two degree's of freedom for which all periodic orbits are 
isolated and unstable. The zeta function then reads 

cxo , —ikT \ '"+1 



p m— 



where Ap is the expanding eigenvalue of Mp. 

Putting w — 1 we obtain a zeta function whose zeros yields the so called resonance spectrum 
or spectrum of correlation exponents. 



_,-i.T \ ™+l 



p m=0 ^ I PI p 
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The leading zero, ko, is the escape rate, which is zero for a bound system. 

By using different weights w one can probe different properties of the flow. 

The topological entropy is obtained by considering the weight uj — |A(a;,i)|. A{x,t) is the 
expanding eigenvalue of the Jacobian transverse to the flow. It is only approximately multiplica- 
tive along the flow but it is possible to modify it slightly so as to become exactly multiplicative 



14 . However, this subtle difference is of no importance for our purposes. The trace formula then 



becomes 

oo 



n=l 



The leading zero lies on the negative imaginary axis. The topological entropy is then h = iko, 
the asymptotic behaviour of the trace is 

oo 

J2T,J2S{t~nT,)^e^' (7) 

p n—l 

so that the number of cycles with periods less than t is ~ e^* /ht, which is the familiar result. 

The semiclassical case is obtained by considering the weight w = y/\A{x, t)\. We now get the 
following zeta function, 

m+l 

Zscik) =1111(1- I (8) 




This zeta function is called the quantum Fredholm determinant [|3|] and is equivalent to the 
Gutzwiller-Voros |l6| zeta function in the semiclassical limit, and have indeed nicer ana- 
lytical features. We have left out the Maslov indices but it is possible to account for them in the 
weight as well. 

2.2 Symmetry decomposition of Z{k) 

Many systems possess finite symmetries. This enables a factorization of the zeta function. This 
was first shown for the semiclassical case in refs. [0 |l^ and for a more general setting in ref. 



19| . The zeta function Z^^k) in is then written as a product Zw{k) — Yl^ Z^^rik) over r, the 



irreducible representations of the symmetry group G, and each Z^, is given by 

oo y _^rp J, , m+l 

Z^Ak) - n n 1 - Xr{9p)Wp^--^ , (9) 

where p runs over all prime cycles in the fundamental domain. Xridp) is the group character for 
the irreducible representation r, and gp is the group element obtained as the product of group 
factors associated with the reflections of the domain walls. 

Throughout the paper we will restrict ourselves to one-dimensional irreducible representations 
of the group. Periodic orbits running along symmetry lines must be given a special treatment. 
We will discuss this when the result is needed. 

In the quantum case the different irreducible representations corresponds to the different 
symmetry classes of the wave functions. But in the general classical case it is harder to give a 
comprehensible interpretation. The important thing to keep in mind in the following is that the 
leading zero is always in the symmetric representation. 

2.3 Zeta functions in the BER approximation 

In ref. an approximate expression for the zeta function is given for intermittent, ergodic 
Hamiltonian systems. The idea is based on a paper by Baladi, Eckmann and Ruelle [ pO| and it 
is therefore referred to as the BER approximation. 
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In an intermittent system there are two, more or less, distinct phases; one regular and one 
irregular (chaotic). Call the consecutive instants when the system enters the regular phase {ti} 
and consider the intervals li = [U-i, ti\. Provided the chaotic phase is chaotic enough, the motions 
in different intervals are nearly mutually independent. This is called assumption A. In particular, 
the lengths of these intervals A; = ti — ti^i are mutually independent and A may be considered as 
a stochastic variable with probability distribution p(A). Under this assumtion, the zeta functions 
may be expressed in terms of the Fourier transform of p{A) 



We use the " - symbol to denote all quantities in the BER approximation. 

This approximation is not restricted to unit weight. If the weight associated with one interval 
is a function of the length of the interval A only, then 



This is one of the central ideas of this paper. To show this, one just repeats the calculation of 
ref. 1^, squeezing in the weight at the appropriate places. As it is straightforward to do we will 
not perform it here. The assumption above, that the weight is a function of A only, is quite 
reasonable for many systems. For the Sinai billiard, which we will study soon, the expression has 
to be slightly generahzed. 



The purpose of this section is to apply the BER approximation, that is formulate Zw,r{k), for 
the Sinai billiard ||2l[| and compare with cycle expansions. The billiard consist of a unit square 
with a scattering disk, having radius i?; < 2R < 1, centered on its midpoint, cf. fig. lb. 
The Sinai billiard is fairly simple but have all the typical features for bound chaotic system; it 
is intermittent and lacks a simple symbolic dynamics. Moreover, it has marginally stable orbits. 
However, a trajectory in the hyperbolic phase (hyperbolic in the sense that all cycles there are 
unstable and isolated) will never end up on a marginally stable orbit, it can only come arbitrary 
close to. This means that we can separate out this marginal phase. That is, the integral (^) is 
only performed over the hyperbolic part of phase space, thus excluding a region of measure zero, 
so that the sum in eq. (|^) runs only over isolated periodic orbit. 

This separation cannot be performed in the quantum case due to an extra contribution from 
the marginal orbits only disappearing in the limit E oo. The study of the quantum case will 
be postponed to a forthcoming paper. 

The trajectory of the Sinai billiard consists of laminar intervals, when bouncing between 
the straight sections, interrupted by scatterings on the central disk. The Sinai billiard seems 
therefore ideally suited for the BER approximation. When the disk radius is small, the memory 
of the previous laminar interval should be almost completely lost. We will not discuss corrections 
to the BER approximation arising from the correlations between laminar intervals. We simply 
assume that they are small when the disk radius is small. The length of the chaotic interval is 
infinitely short so that A is simply the length of the trajectory between two disk scatterings. For 
big radii, there is another source of intermittency [p2| [2^ , the trajectory may be trapped between 
the disk and the straight section. For these reasons we will focuse on the limit of small disk radii. 

3.1 Symbolic dynamics 

First we will define a symbolic dynamics for the Sinai billiard. The reason is twofold. First 
we need it for the application of the BER approximation. Secondly we will use it for finding 
the periodic orbits of the system, which we will need to test the results obtained from the BER 
approximation. 




(10) 




(11) 



3 The Sinai billiard 
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Figure 1: Three equivalent representations of the Sinai billiard, a) The fundamental 
domain, b) The original Sinai billiard with definitions of the variables 4> and a. c) 
The unfolded system with a sample of orbits which are periodic of period one in the 
fundamental domain. 
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We let the disk be our Poincare surface of section. The canonical variables are {£,,p^) — 
{R4>, \/2Esina) where the angles Q < (j) <2tt and a: —j<a<^ are defined in fig. lb. The 
area preserving map i— *■ (CPf) has uniform invariant density. 

We now want to define a symbolic dynamics, such that each iterate of the map corresponds 
to one symbol, that is, each symbol corresponds to one laminar interval. 

To do this we define two equivalent billiard systems, derived from the original system. The 
original billiard has the symmetry group G = C4V, [p4[ . This group has eight elements in five 
conjugacy classes, see table 1. 

The first derived system is the fundamental domain, or the desymmetrized system, see fig la. 
This system has, by definition, no symmetry, except for time reversal. 

To obtain the second derived system we go in the opposite direction. We unfold the billiard, 
see fig Ic, into a regular lattice of disks. Apart from the group C^v this system has also a discrete 
translational symmetry. 

We note that any trajectory in the unfolded billiard may be encoded by a sequence of column 
vectors . . . Qi^iQiQi+i . ■ ■ 

^■{(^ ^ ^ e Z,gcd{n^,ny) = 1} (12) 

where gcd{nx,ny) stands for the greatest common divisor of Ux and Uy. The number [uy) 
simply tells the number of disk sites the segment of the trajectory has travelled in the horisontal 
(vertical) direction. It is easy to see that non comprime Q's cannot be realised. 

A general trajectory corresponds to a semi infinite string of Q's. We will focuse on periodic 
strings Q1Q2 ■ ■ ■ Qn- Such a string uniquely determines the trajectory in the unfolded billiard up 
to a translation which means that the corresponding trajectory is completely determined in the 
original billiard, and is indeed periodic there. Of course, some trajectories defined in this way 
would need to go through disks. The corresponding symbol sequences are then said to be pruned. 
But when the disk radius is small we expect this pruning to apply only to those orbits with some 

long segment(s). More exactly if all the Q's have + <C 1/i? we expect the pruning to be 

neglible. This is very convienient for computing puposes since a large fraction of the generated 
symbol sequences will be realised as true orbits. 

We observe that a symmetry transformation of the orbit corresponds to a transformation 
Qi ^ gQi where g is an element of di^y represented by a 2 x 2 matrix, see table 1. 

Our goal is to encode cycles in the fundamental domain. Due to symmetry, several cycles 
of Q's correspond to one and the same cycle in the fundamental domain. We therefore want 
to translate the cyclic symbol string Q1Q2 ■ ■ ■ Qn to another string S1S2 . . . s„ where the s's are 
taken from another alphabet. The translation must fulfill the following conditions: 

1. The translation is unique. 

2. The sequence S1S2 ■ ■ • s„ is invariant under a symmetry operation on Q1Q2 ■ ■ ■ Qn ■ 

3. A cyclic shift on the Q's corresponds to a cyclic shift on the s's. 

This will be (almost) achieved in the following way. Let each Si be an ordered pair Sj = {qi,gi) 

e w (13) 
w : { ^ ^ ^ ;nx,ny e N,nx > Uy, gcd{nx, Uy) = 1} 

The translation is defined by: 

Qi = .9051 91 

Q2 = 909192 92 (14) 
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Identity 
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Inversion 


R+,R- 


/ -1 \ 




\\ 

I -1 , 




Rotation in the positive and negative directions 


(T| , (T_ 




-1 \ 




r 1 \ 

V -1 ^ 




Reflection in the vertical and horisontal axis 






-1 ^ 

1 ^ 


1 


( \ ^ 

^ -1 , 




Reflections in the diagonals 



Table 1: The elements of the group C4^„. They are represented by 2 x 2 matrices operating 
on the column vectors qi and Qj. 



Qn = 503l ■■■gnq-n 

(15) 

The translation is unique as long as no vector qi is equal to ( J ) or ( |[ ) (we refer to these as 
boundary symbols). We will deal with this problem in a while. Condition 2 is also fulfilled - a 
symmetry operation affects only go. It is easy to check that also condition 3 is fulfilled (a cyclic 
shift also affects go). 

Suppose now that qj is a boundary vector. If we infer the restriction that gj must not be 
a reflection (a\, (t_, a/ or cr\) the translation is again unique provided that not all the g's are 
boundary symbols. 

Suppose the latter is the case and that we have found one translation of the string Q1Q2 ■ ■ ■ Qn, 
namely: S1S2 ■ ■ ■ s„. It is now possible to find another s-string corresponding to the same Q-string. 
This is done in the following way. Suppose that qi has the symmetry (j)i, that is, if qi — { J ) 

then (pi = (J- and if gi = ( \ ) then (pi = a / . Then the transformed sequence 

qi qi 

gi 1-^ (j)i-igi(j)i i ^ (16) 
50 ^ go4>n 



corresponds to the same Q string. Remember that the string is periodic so that e.g. (po = (pn- 
Performing this transformation once more gives back the original s string. So there remains a 
problem that an orbit in the fundamental domain might correspond to two symbol strings and 
we have to account for this problem when generating the periodic orbits. The symbolic dynamics 
defined is still a good one and the problem above should be intepreted as a boundary effect. 

As we said, the problem arises only if all vectors are on the boundary. If one is not, and we 
try to perform the transformation (|l^) , we see that at least one gi corresponding to a boundary 
qi would transform to a reflection and such strings have already been forbidden. 

The group characters in eq. (^ are determined both from the disk bounces and the bounces 
on the square walls. They are easily expressed in therms of the symbol code 

Up 

X.(5p) = nx.(<?.)x.(^-r^-+"«- , (17) 
1=1 

To obtain this we have used the familar property of group characters: Xr{gig2) = Xr{gi)Xr{g2)- 
where gi, nx,i and nx,i are taken from the symbols string corresponding to prime orbit p. 
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One may also consider the Sinai billiard with periodic boundary conditions. One then simply 
exclude the reflection factor Xr{<^-) in cq. 



3.2 Free Directions 

For any disk radius < 2R < 1 there are a finite number of directions through which a trajectory 
may go without ever bouncing on any disk. We call them free directions. Consider the direction 
vector {x,y). A direction cannot be free if x/y is irrational. So we may take {x,y) as coprime 
integers such that x/y < 1. The direction [x, y) is found to be free if 

2R < ^ . (18) 

To show this, consider an integer lattice of points (or consult ref [^). We are seeking the lattice 
point {^,r]) closest to the line from (0,0) to {x,y) and lying below. The closest point above this 
line is then, due to symmetry, {x — ^,y — rj). By the definition of Farey sequences, | is the 
preceeder of - in the Farey sequence of order x. The Farey theorem then says that ^y — rjx = 1 



25 . This is nothing but the area of the parallellogram spanned by the points (0, 0), (^, t]), (x, y) 



and {x — ^,y — rj) so that the distance between the point (^, r]) to the diagonal is 1/ \/ x^ + y"^ and 
the result above follows. 

The direction {x, y) being free is equivalent with the existence of the periodic orbit s where 
s = : ).E). 



3.3 The zeta function in the BER approximation 

We begin by considering the symmetric representation r — Ai (then Xr{gp) = 1) ^-nd unit weight 
It; = 1 in the zeta function. The corresponding indices will be omitted in the zeta function. 

Utilizing the uniformity of the invariant measure pTf we can write down the following expres- 
sion for the probability distribution J3(A) 

p(A)«^a,(0)(5(A-rj , (19) 

where 0^(0) as given by 

ag(0) = - / d(j) d[sin a) . (20) 

is the part of phase space (i.e. the a, 0-plane) for which the trajectory hit disk q (in the 
unfolded system). We have used the symmetry by letting q £ lu, meaning that we take only 
orbits going out into the first octant into account. Only the fractions of the disks ( q ) and ( I ) 
lying in this octant are included. The argument (now put to zero) in ag(0) is introduced for later 
purposes. 

We have made the approximation that all trajectories to disk q have the same length Iq 
corresponding to the time of flight Tg (this will be refined a little in section 3.5). This is reasonable 
when the disk radius is small. Taking the Fourier transform we get: 

Z(fc) = 1-^0,(0)6-*'=^'' . (21) 

Although we use equality sign above we must keep in mind that this result involves an approxi- 
mation in addition to the BER approximation. 



3.4 Relation to cycle expansion of Z{k). 

We are now interested in the relations between the approximate zeta function Z{k) and the cycle 
expansion of the exact zeta function Z{k). First we elaborate a little more on Z{k). 
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Cycle expansion 
BER approximation 



Figure 2: Comparison of the prefactors of eqs (22) and (25) plotted versus sinq 



It is convenient to write ag(0) = Vq^i where ^ is the normalized phase space area taken 
up by disk q if it is not shadowed, and fjq is a pruning factor; < 77^ < I. We thus write 



(22) 



The direction (1,0) is always free if 2R < 1. This means that every disk ( " ) is accessable 
and takes up a normalized phase space area 



a ( .{Q)^-{2R+\/2R-2)- — 

q={ " TT n(n— l)(n — 2) 



(23) 



Let's now consider the zeta function Z(k). First we will neglect the factors with to > in 
(^. Given a periodic orbit S1S2 • ■ ■ Sn we make the following assumtion {assumtion B): 

1. The accessable values for s,; does not depend on Sj, j ^ i. 

2. The stability eigenvalues are given by A ^^...^,^ = Hi ^~ 

3. The periods are given by Tg^ — T— 

This assumption is not at all true, although it is realistic if all Tj's are small. We introduce it as 
a working assumption to shed some light on the nature of the BER approximation. 

Condition 1 means that we can make a curvature expansion Conditions 2 and 3 means 
that all curvatures disappear and what remains is the fundamental part 



(24) 



?7s is a pruning factor which is 1 if cycle s exist, and zero otherwise. Since we have assumed 



that T- 



is a function of q only (= Tq) we sum over g e C'4,y keeping q 



R 



Z(k) - 1 - E ''9^ (1 + \/2)(l + COS(t)q) + sin4>q 



-ikT„ 



) hxed and get 
(25) 
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The angle (j)q is the angle between the horisontal axis and the direction vector {nx,ny). rjq takes 
on discrete values and is = 1 if aU g' s are realised. In fig. 2 we compare the prefactors (in 
arbitrary units) of eqs. (|2^ ) and (|2^) provided that 77, = T)q. The similarity with eq. (^4|) is 
obvious and we may deduce that Assumptions A and B are similar. The reason for the agreement 
is that the available phase space region in eq. (|2^) is extrapolated from the local stabilities in eq. 
(p5|). This should work well as long as there is no shadowing (=obstructing) disks in between. 

We will now see that eq. ( p5|) fails completely when pruning becomes essential, that is when 
fjq <^ 1. This is e. g. the case far out in the free directions. Consider the free direction (1,0). 



We see from (23) that a „ (0) decays as ~ There is only one period-one cycle with 

( 1 ) 



q — { 1 ) , and it is the one having g = and the corresponding prefactors 1 / A-^^-^-^ decays as 
l/r?\ 

The reason for the discrepancy is the fake assumption B. An apparent example of violation 
of Bl is the following. The cycle s where s = {q,g = E) is pruned but the same symbol s may 
very well appear in longer cycles. 

As we have already said, there is no well defined fundamental part of the expanded zeta 
function. Indeed, by our discussion we are led to suggesting that the BER approximation provides 
a natural generalization of the fundamental part, with the extra advantage of preserving unitarity: 
fco = is by construction a zero of the unweighted zeta function Z{k). 

We have seen that the problem of establishing a fundamental part of the zeta function using 
periodic orbits is connected with the marginal orbits which we have pruned by excluding the 
marginal phase and the ackumulation of hyperbolic orbits towards it. The BER approximation 
tells us that there may be simpler means of exploring the available phase space for a certain 
symbol than letting periodic orbits explore it. 

3.5 Other representations, other weigths in the BER approximation 

It is straightforward to study other than the symmetric representation Ai in the BER approxi- 
mation. In this approximation we only consider one single disk bounce. Since we have decided 
that the trajectory goes out into the first octant it suffices to determine what octant it came from 
in order to determine the group element gs associated with this particular bounce. That is, we 
divide the phase space region Q,q into eight regions fis where s — {q,g) and g £ C^^y. Some of 
these may of course be empty. 

The group character associated with the whole interval (one bounce -|- the trajectory to the 
next disk) is then (cf, (0)) 

Xr-(5s)=Xr(5)Xr.(fT-)"-+"'' ■ (26) 

The zeta function may then be written 

Zrik)^l~ J2 X.(ff.)as(0)e-^'=^' , (27) 

S = (9:9) 



where in 0^(0) is defined by ( pO[ ) but the integral extends only over Qs- 

It is also straightforward to use other weigths w in the BER approximation for the Sinai 
billiard. All weights we are going to study in this paper belong to the family w — |A(0, a)\'^ . The 
corresponding family of zeta functions then reads 

ZrAk)-l~Y.^r{9s)as{T)e-'''^^ , (28) 

s 

where 

asir) = 

7^ Jn, 

11 

A(0,a)«— ^ . (30) 
Kcos a 

As A(0, a) is not constant over Qs: c.f. eq. ([ll|), we have averaged it over fig. 



a,{T)^- [ \A{(l),a)\''d<j)d{sina) (29) 
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Figure 3: The trace of the unweighted evolution operator. The parameters are R = 0.1, 
r = Ai,a = 0.1 

4 Computations and results 

We now turn to explicit computations. The aim is to check the validity of the BER approximation 
as outlined in section 3.5, meaning that we restrict ourselves to the family of weights w = \A\'^ . 
The strategy is the following. First, we compute the trace in the BER approximation, 

^^^..^^fe-m^e-^-^/^rffc . (31) 

Then, in order to check this result, we compute the 'exact' trace, by means of the trace 
formula 

^r4..,. - E E X.(gp)"|Ap|- _ _ , , (32) 



To this end we need the periodic orbits. We use gaussian smearing in order to avoid the delta 
peaks. 

The full trace of the full system (fig lb) is obtained as a sum of traces for all irreducible 
representations of G = di.y. As we will not treat two dimensional representations we will not 
perform the sum. 

We mentioned before that orbits running along symmetry lines has to be given a special 
treatment. There are two such cycles in the Sinai billiard. We treat them in the following way. 
One has symmetry cr_. We simply exclude from the trace formula (|3|) if Xr(c_) = —1. The 
other one, having symmetry a/ is treated analogously. This procedure is not exact, c. f. [p^ |l9| , 
but very accurate. 

The computational task lies in determining the quantities as(r) as defined by eq. (p^. 

We do this in the simplest thinkable way by binning the phase space corresponding to a disk 
q and deleting the phase space shadowed by nearer disks. More exactly, we bin the a variable, 
keeping track of the limits in the direction. 

The integrals will for obvious reasons diverge if r > 2. 
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Of course, we only calculate a finite number of as. Then there are infinite tails in all free 
directions. Let us for a minute restrict ourselves to the case r — Ai and the (1,0) free direction. 
Far out in the free directions only a small fraction of disk g = ( " ) is visible. Then ls=(q,g) is 
essentially independent of 5. If we sum over g keeping q fixt 

M 2 2n {2lR-^)^-^'\l-2R) 1 

1 2-T ni+-/2(n-l)i-/2(n-2)i--/2 ^ ' 

We will not bore the reader with details about the other r's and free directions. These infinite 
tails will make the zeta functions diverge in the upper half fc-plane. It is thus crucial that we find 
an analytical continuation of the zeta function. 

We also have to choose the lengths Is associated with $7^. We simply choose it to be the 
period of the cycle s regardless if it is pruned or not. We choose the energy such that the lengths 
Is and the times Ts coincides. 

We will also calculate all periodic orbits up to some maximal period. We thus generate symbol 
codes according to the scheme outlined in the previous section. The orbits are then found by an 
extremum principle according to ref. H^. Last it is checked if each cycle is allowed or pruned. 

4.1 The unweighted zeta function, r = 0. 

We begin by considering the summation of the tails. Consider disk q — { " ) in the free direction 
(1, 0). Only a small fraction of this disk is visible for large n. Then ^s=(q,g) is nonempty only for 
g = E and g = cr^ and all points lie on about the same distance Ig — Iq- Restricting ourselves to 
r = ^1 we have to consider the following series 



-e-'"'" . (34) 

7M //. III 11. /. ' 

n=N+l 



^ ^ n(n — l)(n — 2) 



If we as an approximation put lq^(^n,i) ~ the series reduces to a Fourier series in z = expl—ik). 
This is easily summed since 

nin - 1)(. - 2) - -2^' ' ' + 4^^ - ' ^ ~ + i • (3^) 

The zeta function will thus have branch cuts along Re{k) — 2ttN and Im{k) > 0. 

There are similar contributions in all free direction. The free direction [x, y) will induce 
branch cuts along Re{k) = 2ttN / \J ^ -ip- and Im{k) > 0. 

Close to the origin fc = the zeta function behaves as Z{k) — ci{ik) + d2{ik)'^ \og{ik) + 
C2(ifc)2 . . (r — Ai). Consequently the trace will have the asymptotic behaviour 

i7^4=o..=A,,.^l-CA (36) 

where C is a positive constant, it may be related to the expectation value of p(A) ||]. The 
convergence of the trace formula towards unity is usually reffered to as a periodic orbit sum rule 
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When considering other representations, r, some tail series will have alternating signs. This 
will shift the branch cuts in the real direction. These zeta functions will not have a root at the 
origin and the traces goe asymptotically to zero. 

We now turn to the numerical computation of the trace. When we calculate the Fourier 
transform we simply make the integration along the the real axis, except at the origin which we 
sneak just below. 

With the number of periodic orbits calculated (between 5000 and 8000 depending on R) 
we cannot check the asymptotic limit discussed above. To do the numerical check we restrict 
ourselves to much shorter times t. The result is seen in figs 3 and 4. The agreement is in most 
cases very good. We note 
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Figure 4: Trace of the unweighted evolution operator obtained from the BER approxi- 
mation and directly from periodic orbits 
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• For small t the agreement is less good or even lousy. The exact result for small t depends 
on a small number of cycles and the BER theory requires an average effect from many 
cycles. 

• For fixed representation r, the agreement is better for small radii R. This is exactly what 
we expected. 

• For fixed radius R the agreement is best for the symmetric representation Ai and worst 
for the anti symmetric ^2- In the former case all Qs provides the same sign but otherwise 
different phase space regions come with different signs (provided by the group characters 
x)- It takes more periodic orbits in order to provide averages for this finer partition. 

We have seen that the trace is almost entirely governed by the branch cuts. (There are a 
few isolated zeros but they are extremely sparse and have no practical importance). We are thus 
faced with an entirely different situation than for Axiom-A systems for which the zeta function 
is entire and the zeros gives the discrete spectrum of the operator. 

We have not discussed in what sense the cuts may be interpreted as a continous part of the 
spectra. If the trace, by means of the sum over eigenvalues, is not mathematically well defined, 
we let it be defined by means of the corresponding periodic orbit sum. 

4.2 Topological entropy, r = 1 

We now have to find the analytic continuation for sums of the form, cf. cq. (^3|), 

-ikcn 



S ^3/2 
n=N+l 



(37) 



This is just a power series in exp{~ikc) . We observe that 



n— 1 ^ ^ 



Using Stirlings formula we find the following asymptotic expression for the coefficient 

(2n-3)!! 1 1 
(2n)!! ^ V4^n3/2 



(39) 



enabling an summation of the (|37|), provided of course that N is sufficiently big. The zeta function 
still have cuts in the upper half plane. However, there are isolated zeros far down in the lower 
half plane. The determination of these does not depend very much on the summation of the tails. 
Therefore we could take only the leading order of the tails series into account above. 



When preforming the integral (31) we let the contour go below all singularities of the inte- 



grand, that is below all discrete zeros ka of Z. The sum over residues gives 

= l,r=Ai,cr ~ ^ 



trCU^^.^^ye^'^''^e~'^>'/' . (40) 



The result is seen in fig. 5. This result is again good for large t. 

Of course, in the contour integration we have neglected the contribution coming from the 
integration around the branch cuts cf ref. This should give an extra contribution ~ — l/\/t 
which is of course overwhelmed by the exponential increase, but could explain the discrepancy 
for small t in fig. 5. 

We have determined the topological entropy h{R) for two different radii R. We find h{0.2) — 
2.1346 and h{0.1) = 2.02024. The entropy h{R) tends to a finite limit when R ^ 0. It hmiting 
value is obtained from the equation 

1- e-''(°)V^ = o . (41) 

i,jeZ, gcd{i,j) = l 
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Figure 5: The trace of the topological evolution operator obtained from the BER ap- 
proximation and directly from periodic orbits 

This is derived directly from eqs. ( p8|) and (^9|) but is indeed exact, thus demonstrating consis- 
tency of our approach. It allows determination of h{0) within a fraction of a second of computer 



time. W e find h{0) = 1.9133307629.... This value differs slightly from Berry's estimate [H 
h{0) « ^127^= 1.95441. 



Considering the one parameter family of weights |A|'^, we see that if r > the leading zeros 
lie below the branch cut and the trace exhibit exponential behaviour. When r = the leading 
zeros coincides with a branch point. When t < the behaviour is entirely given by the branch 
cuts, i.e. it is goverened by the tails in the free directions. This is usually referred to as a phase 
transition We have thus been able to give a nice analytical description of such a phase 

transition in terms of the zeta function; as the collision of the leading zero with a singularity. 

5 Discussion 

The aim has been to show that the BER approximation works, and works well for quite general 
weigths w. We have not aimed at calculating more physically interesting properties, such as decay 
of correlations, fractal dimensions etc. This is of course a natural next step. 

Regarding correlation functions one may, loosely, say that they behave very much like the 
the trace of the evolution operator with unit weight. This is easily seen if the eigenvalues 
exp(ifc„t) are isolated. Then the correlation function C{t) may (formally) be expanded as C{t) = 

c„ exp(iA:„<). The coefficients c„ depend on the observables. The trace trC* — ^^exp(iA:„t) 
may therefore be thought of as a archetype correlation function. The velocity autocorrelation 
function for the Sinai billard (with continous time) have been reported to decay as 1 /t , which 
is the same decay as we found for the trace. 

The BER approximation is in principle applicable to a wide range of bound ergodic Hamilto- 
nian systems. The details have to be worked out for each separate case. For the Sinai billiard it 
turned out to be particularly simple. The reason for this is that we could make the approximation 
that the lengths Is are constant over the phase space regions ilg- This was because the flow is 
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suspended in the sense that hminar segments of arbitrary lengths does not exist. For instance, 
there is no such length shorter than 1 — 2R. This has to do with the presence of marginally stable 
orbits. A consequence of this is the periodicity of the cuts. 

The situation is very different for e.g. the hyperbola billiard. In ref. 1^ we found evidence 
that the zeta function (r — 0) for this system has a discrete spectrum of zeros and a single 
cut along the positive imaginary axis. We used cycle expansions to obtain this result, but since 
the zeta function is divergent we had to do some ad hoc manipulations. The expansion used 
exhibited some obvious similarities with the BER approximation worked out in this paper. The 
expansion consisted almost exclusively of cycles with only one laminar segment. To avoid (the 
most serious) divergence we selected out an infinite subsequence ackumulating towards a pruned 
orbit and responsible for the branch cut. This sequence has as a counterpart the tails in the free 
directions of the Sinai billiard. 

Perhaps the most powerful result in this article is the possibility of making asymptotic state- 
ments, cf eq (^), about periodic orbits. By varying the parameter r one can obtain a wealth of 
such periodic orbit sum rules. This asymptotics is essential for doing semiclassical calculation on 
spectral fluctuations based on Gutzwiller's trace formula ^. We have therfore outlined how 
to do periodic orbit theory without a single periodic orbit. 
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